Augmenting Biogas Process Modeling by Resolving Intracellular Metabolic Activity

The process of anaerobic digestion in which waste biomass is transformed to methane by complex microbial communities has been modeled for more than 16 years by parametric gray box approaches that simplify process biology and do not resolve intracellular microbial activity. Information on such activity, however, has become available in unprecedented detail by recent experimental advances in metatranscriptomics and metaproteomics. The inclusion of such data could lead to more powerful process models of anaerobic digestion that more faithfully represent the activity of microbial communities. We augmented the Anaerobic Digestion Model No. 1 (ADM1) as the standard kinetic model of anaerobic digestion by coupling it to Flux-Balance-Analysis (FBA) models of methanogenic species. Steady-state results of coupled models are comparable to standard ADM1 simulations if the energy demand for non-growth associated maintenance (NGAM) is chosen adequately. When changing a constant feed of maize silage from continuous to pulsed feeding, the final average methane production remains very similar for both standard and coupled models, while both the initial response of the methanogenic population at the onset of pulsed feeding as well as its dynamics between pulses deviates considerably. In contrast to ADM1, the coupled models deliver predictions of up to 1,000s of intracellular metabolic fluxes per species, describing intracellular metabolic pathway activity in much higher detail. Furthermore, yield coefficients which need to be specified in ADM1 are no longer required as they are implicitly encoded in the topology of the species’ metabolic network. We show the feasibility of augmenting ADM1, an ordinary differential equation-based model for simulating biogas production, by FBA models implementing individual steps of anaerobic digestion. While cellular maintenance is introduced as a new parameter, the total number of parameters is reduced as yield coefficients no longer need to be specified. The coupled models provide detailed predictions on intracellular activity of microbial species which are compatible with experimental data on enzyme synthesis activity or abundance as obtained by metatranscriptomics or metaproteomics. By providing predictions of intracellular fluxes of individual community members, the presented approach advances the simulation of microbial community driven processes and provides a direct link to validation by state-of-the-art experimental techniques

Genome and proteome analyses show the gaseous alkane degrader Desulfosarcina sp. strain BuS5 as an extreme metabolic specialist

© The Author(s), 2022. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Chen, S.-C., Ji, J., Popp, D., Jaekel, U., Richnow, H.-H., Sievert, S. M., & Musat, F. Genome and proteome analyses show the gaseous alkane degrader Desulfosarcina sp. strain BuS5 as an extreme metabolic specialist. Environmental Microbiology, 24, (2022): 1964-1976, https://doi.org/10.1111/1462-2920.15956.The metabolic potential of the sulfate-reducing bacterium Desulfosarcina sp. strain BuS5, currently the only pure culture able to oxidize the volatile alkanes propane and butane without oxygen, was investigated via genomics, proteomics and physiology assays. Complete genome sequencing revealed that strain BuS5 encodes a single alkyl-succinate synthase, an enzyme which apparently initiates oxidation of both propane and butane. The formed alkyl-succinates are oxidized to CO2 via beta oxidation and the oxidative Wood–Ljungdahl pathways as shown by proteogenomics analyses. Strain BuS5 conserves energy via the canonical sulfate reduction pathway and electron bifurcation. An ability to utilize long-chain fatty acids, mannose and oligopeptides, suggested by automated annotation pipelines, was not supported by physiology assays and in-depth analyses of the corresponding genetic systems. Consistently, comparative genomics revealed a streamlined BuS5 genome with a remarkable paucity of catabolic modules. These results establish strain BuS5 as an exceptional metabolic specialist, able to grow only with propane and butane, for which we propose the name Desulfosarcina aeriophaga BuS5. This highly restrictive lifestyle, most likely the result of habitat-driven evolutionary gene loss, may provide D. aeriophaga BuS5 a competitive edge in sediments impacted by natural gas seeps.This study was financed by the Max Planck Society and by the Helmholtz Association of German Research Centres. The draft genome was sequenced as part of the U.S. Department of Energy Joint Genome Institute (DOE-JGI) Community Science Program project 1078203 awarded to S. M. Sievert and F. Musat. The work conducted by the DOE-JGI, a DOE Office of Science User Facility, is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231. Lynne Goodwin (Los Alamos National Laboratory) is acknowledged for project management support of the draft genome sequencing. Further support was provided by the U.S. National Science Foundation grant MCB-0702677 (to SMS), and by the Helmholtz Association grant ERC-RA-0020 (to FM). We acknowledge the Centre for Chemical Microscopy (ProVIS) platform at the Helmholtz Centre for Environmental Research – UFZ, for using their analytical facilities. ProVIS is supported by European Regional Development Funds (EFRE – Europe funds Saxony). We acknowledge the Bundesministerium für Bildung und Forschung (BMBF)-funded German Network for Bioinformatics Infrastructure de.NBI (031A537B, 031A533A, 031A538A, 031A533B, 031A535A, 031A537C, 031A534A, 031A532B) for providing computational resources

Ammonia Inhibition of Anaerobic Volatile Fatty Acid Degrading Microbial Communities

Ammonia inhibition is an important reason for reactor failures and economic losses in anaerobic digestion. Its impact on acetic acid degradation is well-studied, while its effect on propionic and butyric acid degradation has received little attention and is consequently not considered in the Anaerobic Digestion Model No. 1 (ADM1). To compare ammonia inhibition of the degradation of these three volatile fatty acids (VFAs), we fed a mixture of them as sole carbon source to three continuous stirred tank reactors (CSTRs) and increased ammonium bicarbonate concentrations in the influent from 52 to 277 mM. The use of this synthetic substrate allowed for the determination of degradation efficiencies for the individual acids. While butyric acid degradation was hardly affected by the increase of ammonia concentration, propionic acid degradation turned out to be even more inhibited than acetic acid degradation with degradation efficiencies dropping to 31 and 65% for propionic and acetic acid, respectively. The inhibited reactors acclimatized and approximated pre-disturbance degradation efficiencies toward the end of the experiment, which was accompanied by strong microbial community shifts, as observed by amplicon sequencing of 16S rRNA genes and terminal restriction fragment length polymorphism (T-RFLP) of mcrA genes. The acetoclastic methanogen Methanosaeta was completely replaced by Methanosarcina. The propionic acid degrading genus Syntrophobacter was replaced by yet unknown propionic acid degraders. The butyric acid degrading genus Syntrophomonas and hydrogenotrophic Methanomicrobiaceae were hardly affected. We hypothesized that the ammonia sensitivity of the initially dominating taxa Methanosaeta and Syntrophobacter led to a stronger inhibition of the acetic and propionic acid degradation compared to butyric acid degradation and hydrogenotrophic methanogenesis, which were facilitated by the ammonia tolerant taxa Syntrophomonas and Methanomicrobiaceae. We implemented this hypothesis into a multi-taxa extension of ADM1, which was able to simulate the dynamics of both microbial community composition and VFA concentration in the experiment. It is thus plausible that the effect of ammonia on VFA degradation strongly depends on the ammonia sensitivity of the dominating taxa, for syntrophic propionate degraders as much as for acetoclastic methanogens

Determination of Microbial Maintenance in Acetogenesis and Methanogenesis by Experimental and Modeling Techniques

For biogas-producing continuous stirred tank reactors, an increase in dilution rate increases the methane production rate as long as substrate input can be converted fully. However, higher dilution rates necessitate higher specific microbial growth rates, which are assumed to have a strong impact on the apparent microbial biomass yield due to cellular maintenance. To test this, we operated two reactors at 37°C in parallel at dilution rates of 0.18 and 0.07 days-1 (hydraulic retention times of 5.5 and 14 days, doubling times of 3.9 and 9.9 days in steady state) with identical inoculum and a mixture of volatile fatty acids as sole carbon sources. We evaluated the performance of the Anaerobic Digestion Model No. 1 (ADM1), a thermodynamic black box approach (TBA), and dynamic flux balance analysis (dFBA), to describe the experimental observations. All models overestimated the impact of dilution rate on the apparent microbial biomass yield when using default parameter values. Based on our analysis, a maintenance coefficient value below 0.2 kJ per carbon mole of microbial biomass per hour should be used for the TBA, corresponding to 0.12 mmol ATP per gram dry weight per hour for dFBA, which strongly deviates from the value of 9.8 kJ Cmol h-1 that has been suggested to apply to all anaerobic microorganisms at 37°C. We hypothesized that a decrease in dilution rate might select taxa with minimized maintenance expenditure. However, no major differences in the dominating taxa between the reactors were observed based on amplicon sequencing of 16S rRNA genes and terminal restriction fragment length polymorphism analysis of mcrA genes. Surprisingly, Methanosaeta dominated over Methanosarcina even at a dilution rate of 0.18 days-1, which contradicts previous model expectations. Furthermore, only 23–49% of the bacterial reads could be assigned to known syntrophic fatty acid oxidizers, indicating that unknown members of this functional group remain to be discovered. In conclusion, microbial maintenance was found to be much lower for acetogenesis and methanogenesis than previously assumed, likely due to the exceptionally low growth rates in anaerobic digestion. This finding might also be relevant for other microbial systems operating at similarly low growth rates

Recently evolved combination of unique sulfatase and amidase genes enables bacterial degradation of the wastewater micropollutant acesulfame worldwide

Xenobiotics often challenge the principle of microbial infallibility. One example is acesulfame introduced in the 1980s as zero-calorie sweetener, which was recalcitrant in wastewater treatment plants until the early 2010s. Then, efficient removal has been reported with increasing frequency. By studying acesulfame metabolism in alphaproteobacterial degraders of the genera Bosea and Chelatococcus, we experimentally confirmed the previously postulated route of two subsequent hydrolysis steps via acetoacetamide-N-sulfonate (ANSA) to acetoacetate and sulfamate. Genome comparison of wildtype Bosea sp. 100-5 and an acesulfame degradation-defective mutant revealed the involvement of two plasmid-borne gene clusters. The acesulfame-hydrolyzing sulfatase is strictly manganese-dependent and belongs to the metallo beta-lactamase family. In all degraders analyzed, it is encoded on a highly conserved gene cluster embedded in a composite transposon. The ANSA amidase, on the other hand, is an amidase signature domain enzyme encoded in another gene cluster showing variable length among degrading strains. Transposition of the sulfatase gene cluster between chromosome and plasmid explains how the two catabolic gene clusters recently combined for the degradation of acesulfame. Searching available genomes and metagenomes for the two hydrolases and associated genes indicates that the acesulfame plasmid evolved and spread worldwide in short time. While the sulfatase is unprecedented and unique for acesulfame degraders, the amidase occurs in different genetic environments and likely evolved for the degradation of other substrates. Evolution of the acesulfame degradation pathway might have been supported by the presence of structurally related natural and anthropogenic compounds, such as aminoacyl sulfamate ribonucleotide or sulfonamide antibiotics

PHIP-associated Chung-Jansen syndrome: Report of 23 new individuals

In 2016 and 2018, Chung, Jansen and others described a new syndrome caused by haploinsufficiency of PHIP (pleckstrin homology domain interacting protein, OMIM *612,870) and mainly characterized by developmental delay (DD), learning difficulties/intellectual disability (ID), behavioral abnormalities, facial dysmorphism and obesity (CHUJANS, OMIM #617991). So far, PHIP alterations appear to be a rare cause of DD/ID. “Omics” technologies such as exome sequencing or array analyses have led to the identification of distinct types of alterations of PHIP, including, truncating variants, missense substitutions, splice variants and large deletions encompassing portions of the gene or the entire gene as well as adjacent genomic regions. We collected clinical and genetic data of 23 individuals with PHIP-associated Chung-Jansen syndrome (CHUJANS) from all over Europe. Follow-up investigations (e.g. Sanger sequencing, qPCR or Fluorescence-in-situ-Hybridization) and segregation analysis showed either de novo occurrence or inheritance from an also (mildly) affected parent. In accordance with previously described patients, almost all individuals reported here show developmental delay (22/23), learning disability or ID (22/23), behavioral abnormalities (20/23), weight problems (13/23) and characteristic craniofacial features (i.e. large ears/earlobes, prominent eyebrows, anteverted nares and long philtrum (23/23)). To further investigate the facial gestalt of individuals with CHUJANS, we performed facial analysis using the GestaltMatcher approach. By this, we could establish that PHIP patients are indistinguishable based on the type of PHIP alteration (e.g. missense, loss-of-function, splice site) but show a significant difference to the average face of healthy individuals as well as to individuals with Prader-Willi syndrome (PWS, OMIM #176270) or with a CUL4B-alteration (Intellectual developmental disorder, X-linked, syndromic, Cabezas type, OMIM #300354). Our findings expand the mutational and clinical spectrum of CHUJANS. We discuss the molecular and clinical features in comparison to the published individuals. The fact that some variants were inherited from a mildly affected parent further illustrates the variability of the associated phenotype and outlines the importance of a thorough clinical evaluation combined with genetic analyses for accurate diagnosis and counselling

A Modular Metagenomics Pipeline Allowing for the Inclusion of Prior Knowledge Using the Example of Anaerobic Digestion

Metagenomics analysis revealing the composition and functional repertoire of complex microbial communities typically relies on large amounts of sequence data. Numerous analysis strategies and computational tools are available for their analysis. Fully integrated automated analysis pipelines such as MG-RAST or MEGAN6 are user-friendly but not designed for integrating specific knowledge on the biological system under study. In order to facilitate the consideration of such knowledge, we introduce a modular, adaptable analysis pipeline combining existing tools. We applied the novel pipeline to simulated mock data sets focusing on anaerobic digestion microbiomes and compare results to those obtained with established automated analysis pipelines. We find that the analysis strategy and choice of tools and parameters have a strong effect on the inferred taxonomic community composition, but not on the inferred functional profile. By including prior knowledge, computational costs can be decreased while improving result accuracy. While automated off-the-shelf analysis pipelines are easy to apply and require no knowledge on the microbial system under study, custom-made pipelines require more preparation time and bioinformatics expertise. This extra effort is minimized by our modular, flexible, custom-made pipeline, which can be adapted to different scenarios and can take available knowledge on the microbial system under study into account

Metabolic Adaptation of Methanogens in Anaerobic Digesters Upon Trace Element Limitation

Anaerobic digestion (AD) is a complex multi-stage process relying on the activity of highly diverse microbial communities including hydrolytic, acidogenic and syntrophic acetogenic bacteria as well as methanogenic archaea. The lower diversity of methanogenic archaea compared to the bacterial groups involved in AD and the corresponding lack of functional redundancy cause a stronger susceptibility of methanogenesis to unfavorable process conditions such as trace element (TE) deprivation, thus controlling the stability of the overall process. Here, we investigated the effects of a slowly increasing TE deficit on the methanogenic community function in a semi-continuous biogas process. The aim of the study was to understand how methanogens in digester communities cope with TE limitation and sustain their growth and metabolic activity. Two lab-scale biogas reactors fed with distillers grains and supplemented with TEs were operated in parallel for 76 weeks before one of the reactors was subjected to TE deprivation, leading to a decline of cobalt and molybdenum concentrations from 0.9 to 0.2 mg/L, nickel concentrations from 2.9 to 0.8 mg/L, manganese concentrations from 38 to 18 mg/L, and tungsten concentrations from 1.4 to 0.2 mg/L. Amplicon sequencing of mcrA genes revealed Methanosarcina (72%) and Methanoculleus (23%) as the predominant methanogens in the undisturbed reactors. With increasing TE limitation, the relative abundance of Methanosarcina dropped to 67% and a slight decrease of acetoclastic methanogenic activity was observed in batch tests with 13C-methyl-labeled acetate, suggesting a shift toward syntrophic acetate oxidation coupled to hydrogenotrophic methanogenesis. Metaproteome analysis revealed abundance shifts of the enzymes involved in methanogenic pathways. Proteins involved in methylotrophic and acetoclastic methanogenesis decreased in abundance while formylmethanofuran dehydrogenase from Methanosarcinaceae increased, confirming our hypothesis of a shift from acetoclastic to hydrogenotrophic methanogenesis by Methanosarcina. Both Methanosarcina and Methanoculleus increased the abundance of N5-methyltetrahydromethanopterin-coenzyme M methyltransferase and methyl-coenzyme M reductase. However, these efforts to preserve the ion motive force for energy conservation were seemingly more successful in Methanoculleus. We conclude that both methanogenic genera use different strategies to stabilize their energy balance under TE limitation. Methanosarcina switched from TE expensive pathways (methylotrophic and acetoclastic methanogenesis) to hydrogenotrophic methanogenesis. Methanoculleus showed a higher robustness and was favored over the more fastidious Methanosarcina, thus stabilizing reactor performance under TE limitation

Approach to Cohort-Wide Re-Analysis of Exome Data in 1000 Individuals with Neurodevelopmental Disorders

The re-analysis of nondiagnostic exome sequencing (ES) has the potential to increase diagnostic yields in individuals with rare diseases, but its implementation in the daily routines of laboratories is limited due to restricted capacities. Here, we describe a systematic approach to re-analyse the ES data of a cohort consisting of 1040 diagnostic and nondiagnostic samples. We applied a strict filter cascade to reveal the most promising single-nucleotide variants (SNVs) of the whole cohort, which led to an average of 0.77 variants per individual that had to be manually evaluated. This variant set revealed seven novel diagnoses (0.8% of all nondiagnostic cases) and two secondary findings. Thirteen additional variants were identified by a scientific approach prior to this re-analysis and were also present in this variant set. This resulted in a total increase in the diagnostic yield of 2.3%. The filter cascade was optimised during the course of the study and finally resulted in sensitivity of 85%. After applying the filter cascade, our re-analysis took 20 h and enabled a workflow that can be used repeatedly. This work is intended to provide a practical recommendation for other laboratories wishing to introduce a resource-efficient re-analysis strategy into their clinical routine